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We consider an inflationary curvaton scenario, where the curvaton decays into two non-interacting 
relativistic fluids and later during the cosmological evolution one of them becomes non-relativistic, 
forming dark matter component of the universe. We study the thermic properties and the generation 
of non-gaussianity in this three fluid curvaton model. By solving the evolution of the system and us- 
ing several cosmological conditions we find that the allowed parameter space is strongly constrained. 
The naturalness of this curvaton scenario is also discussed. 



I. INTRODUCTION 

The most recent data gathered by WMAP [TJ |2j [3] is consistent with the hypothesis that the perturbations in 
the cosmic microwave background radiation (CMB) were generated in an era of cosmic inflation. This is generally 
achieved with a slowly rolling scalar field which leads to an exponential expansion of the universe and the observed 
anisotropy is generated by the fluctuations of this inflaton field jH |5] . Such a minimal scenario leads to adiabatic and 
gaussian perturbations in accordance with the current data. 

A well motivated alternative to the simplest inflationary scenario is the curvaton mechanism [BJ |5J EH EH EJ H2J 
HH HSl |T7l KH HSl 1201 HH [231 [SI 1^1 HH1 HS1 EU1 EH 1321 1^1 EH ESI EE] in which the perturbations 

are generated by a second scalar field, called the curvaton, which stays subdominant during inflation and the actual 
expansion of space is still driven by the inflaton. This allows the inflation potential to have more natural properties 
compared to the single field scenario [T7]- However, the extra degrees of freedom in the system now allow for the 
possibility that the final state is not necessarily purely adiabatic. Instead, the generation of an observable amount 
of isocurvature perturbations is a possibility that can distinguish the curvaton scenario from the simple single field 
inflationary model. 

Early models of curvaton cosmology usually assumed that the decay of the curvaton field produces only radiation 
[Ml ES]- Later models included the possibility of radiation and matter production (SSJ EH]- In the present paper 
we have expanded this line of thought to include the possibility that the curvaton decays into two non-interacting 
relativistic components. This makes it possible to model a cold dark matter component that is initially relativistic 
but turns non-relativistic as the universe expands and dilutes the energy densities. This process can lead to some 
interesting consequences: for example the curvature perturbations of matter are no longer conserved after curvaton 
decay if the decay occurs when the matter component is still relativistic. 

In addition to the background evolution we will also solve numerically the first and second order perturbation 
equations with different initial values and search for physically plausible initial values in a parameter space based on 
the initial values of these variables. Besides studying the correctness of the thermodynamical evolution of the system 
we will also calculate the generated level of non-gaussianity and compare it to the limits given by WMAP 5 year 
data [5] . This strategy makes it possible to study the naturalness of different dark matter particle species within the 
curvaton decay model. 

This paper is organized as follows. In section II we present the relevant quantities and equations of motion of the 
perturbation up to second order. In section III we present the equations of motion of the curvaton model and derive 
the equations of motion for the temperatures for radiation and dark matter. In section IV we present numerical 
results of the simulations and we end this article with discussion and conclusions in section V. 



II. PERTURBATIONS AT FIRST AND SECOND ORDER 

The theory of second order cosmological perturbations has been studied rigorously in the recent years and a good 
introduction to the subject can be found e.g. in |40j . We adopt their notations and conventions, and employ a 
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spatially flat Friedmann-Robertson-Walker -background. The metric tensor up to second order is [10] 



g^dx^dx u = - (1 + 20^ + 4>^)dt 2 + a(i)(^ (1) + -u}f ] )dtdx 



(1) 



dx l dx :> , 



where a(i) is the scale factor and 4>^ r \ Cj^ r \ ip^ and x^ij are (scalar, vector and tensor) perturbation functions at 
first (r = 1) and second order (r = 2) . Written in this form, different gauges can be straightforwardly given in terms 
of the perturbation functions: for example the Poisson gauge is defined as uj^ — x^ij = Xij — an d the spatially 
flat one is V (r) = X (r) = 

A useful set of equations can be derived from the continuity equations = Q\, where T? v is the energy- 
momentum tensor, Qi describes the energy transfer between different fluids and semi-colon subscript (;) denotes the 
covariant derivative. From the continuity equations it follows that the equations determining the background evolution 
of individual fluids are The background evolution equations for individual fluids follow from the continuity equations 

fH = -3H(l + Wi)pi + Qi, (2) 

where LOi — Pi/ pi is the equation of state of th ith fluid and '= d/dt i.e. derivative with respect to physical time. 
The evolution equations of the first order perturbed energy and pressure densities (on large scales) are @T] 

8pf ] + 3H(S P ^ + 5PP>) - 3( Pi + Pi)ip iX) = + SQ^. (3) 

and at second order 

Spf ] + ZH{8pf ] + 5Pl 2) ) - 3( Pi + P^ - 6^[8 Pi + 8Pi + 2( Pi + P t )^] 
= Q l ^ + 8Q^-QM (1) ? + W {1) 5Q<t ) . 

These equations can be simplified further by selecting a gauge. We will be using the spatially flat gauge and the 
Einstein equations simplify [3D] in this case into conditions ip^ — = and 2</>W = —Sp/po at first order and at 
second order the — 0-component of the Einstein equations gives 

(2) = -^ + 4(0«) 2 . (5) 
2 po 



A. Curvature perturbations 

One way to study the evolution of perturbations is to use gauge-invariant curvature perturbations, which relate to 
curvature perturbations on homogeneous-density surfaces. At first order they are defined for component i as 

Cj x) = - (6) 

Pi 

where ' = d/d(1n(a)). At second order the corresponding quantity is defined as 

r (2) r (1)/ r CO r (0 s (0 9 // 

<f } = + 2^V^ + 2^VW W ' + 2^W) - ( S 4^) (^ 2). (7) 

Pi p'i p'i Pi W- / V< > 

Note that we are here neglecting gradient terms since we are only interested in the large scale behavior of perturbations. 

The equations of motion of the first order curvature perturbations can be derived from eq. ^ and Einstein 
equations and this results in 

.(i), ^F mt (,) 8Q int{i) H' Qi 

C > H^-Jl^-^ (8) 

where 8P int(i) = 8P t -p'^pjp^ and 8Q int{l) = 8Q % - Q'fipi/p^. 
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At second order the corresponding equations read as [35] 



- 2 ^ + 4 2 -3 (1+ ^cf>c--4(^ + ^)c'1'-[(^-^^)(cf0l' 

There are instances when the definition of different curvature perturbations might fail e.g. when p\ = 0. Therefore 
in our numerical evaluations we have used the spatially flat gauge and evaluated the density perturbations of different 
components in this gauge at first and second-order. The corresponding equations of motion can be easily read from 
equations (|3| and W) by going to the spatially flat gauge i])^ = % — 0. 



III. THE CURVATON MODEL 



The curvaton model of the present paper consists of three fluids: curvaton, radiation, and a matter component 
that is initially relativistic but becomes non-relativistic with the expansion of the universe. This happens roughly 
when the average energy of matter particles drops below its mass. This non-relativistic matter component can thus 
be thought to form the dark matter of the universe. 

We denote the curvaton by subscript a, radiation by 7 and matter by m. The matter and radiation components 
are assumed to be non-interacting with each other which means there are no processes that would thermalize these 
two components after the curvaton decay. This implies that the curvaton mass is higher than the Hubble rate as it 
starts to decay. Otherwise there would be curvaton mediated thermalization processes that would make the matter 
and radiation components interacting. 

In the current scenario before the curvaton decay the energy density of the universe is dominated by radiation with 
a small component of coherent curvaton field. After the curvaton decay the system consists of two fluids, radiation and 
matter, with separate thermal equilibrium distributions and temperatures T 7 and T mi for the radiation and matter 
component respectively. This is a practical assumption which requires that an unknown self-interactions maintains 
kinetic equilibrium of the matter component. Without this assumption the matter particle distribution should be 
calculated directly from the Boltzmann equations. This would however be quite a demanding task requiring that the 
coherent oscillations of curvaton field are carefully taken into account. Thus we simplify the analysis with equilibrium 
distributions. It should be noted that if there were processes maintaining equilibrium between radiation and matter 
we would enter the WIMP scenario [39 or we would have m ~ T eq < 1 eV and the matter component would begin 
to dominate cosmic evolution immediately after becoming non-relativistic. This is hardly a plausible cosmological 
scenario for the dark matter [47] . 



A. Evolution equations 

For the radiation we have o> 7 = 1/3 but the equation of state of the dark matter component depends on the 
temperature of the matter component, T m . Initially the matter is relativistic and u> m = 1/3 but as the universe cools 
down the matter turns into non-relativistic and uj m — » 0. Because of this T m has time dependence that will be derived 
in the next subsection. This equation naturally needs to be incorporated into the equations of motion presented in 
this subsection. The interaction terms are |38j 

Qcr ^fPa TmPcr 

Q y = TjPa (10) 

where I\ denote the strengths of the interactions. 

Background equations (pi can be written in terms of fractional densities Hi = Pi/p for which the equations of 
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motion are: 



Q! a = Q CT (fi 7 + 3w m fi m ) + 

o., = r^ 7 (ri 7 -t- 3(^ m r^ r 



ffp' 



,+ Hp' 



Q' = 

m 

1 



Qm 



(11) 



l+3« 7 



From the definition of fl; it can be easily seen that Q a + f2 7 + Sl r 
redundant. 

From eq. (|3| we can read the equations of motion for the first order density perturbations 



1, which means that one of the fli equations is 



<5p« , 6QV 



(i) 



if 2p 

7 1 



H 



H 2p 



H 



(12) 



*P«' = -3(l + c, m «>-%^ 



if 



in the spatially flat gauge. At second order the corresponding equations in the flat gauge are 



2Hp 



(* P w) ! 



2Hp 
-3(l + w m )Jp« 



5p< 2 



Fp 
SpW^Q 



2iip 



<5p< 



2) 



ffp ' 



(13) 



H 



Hp 



We also need the gauge invariant curvature perturbations when calculating the non-gaussianity parameter f 1 



NL- 



In the spatially flat gauge these are 



c 



(1) 



(i) 



.(2) 



P'i ' 
(2) 



<5p 



3(1 + ^)](C (1) ) ; 



Qi 



5Q\ 



Hp\ Hp[ 



(14) 



r Qi i Qipt 



PrH 2 p^ffp 



The set of equations (111, ( 12 ) and ( 13 ) can now be evaluated numerically once the initial values have been set. We 
have chosen the system to be initially radiation dominated and non-adiabatic at first and second-order. The initial 
values for the three fluid curvaton decay were derived in |44j and they are given by 



C (1) =1 

C (2) ~l/2 



(15) 



and the latter approximation is valid when T a -C H . The other perturbations are initially set to be zero i.e. 

AT-) = = a(2) 

^7, in ^m,in ^7, in 



Win = Cm in = CyjL = Cm.L = ® since we have assumed that only the curvaton field has initial perturbations. 



4 



B. Thermodynamics of the system 



Our curvaton models starts in a radiation dominated universe with initial radiation and dark matter temperature 
set to T 7i ; n and T m .j n respectively. We will assume that the dark matter particles have a mass m m and by setting 
the initial matter temperature T m ^ n we can change the initial equation of state u> m of the dark matter. Different 
temperatures for radiation and dark matter mean that there are no interactions between them and the temperature 
T m is due to kinetic equilibrium of the dark matter particles. 

We will start by deriving the equations of motion for the variables that govern the evolution of the system i.e. the 



temperatures /3 7 and fi m where fi = T 



The energy density of dark matter particles can be naturally written as 
dp 3 



Ef m = C / d PP 'Ef m ([3 m (E - im)) 



(16) 



where we have made the ansatz /„, = f m (/3 m (E — /i)) for the distribution function and used the notation C = l/(27r 2 ). 
The pressure and number density are defined similarly as 



dp 3 p 2 
(2tt) 3 3E 



fm(fim(E - At)) 



and 



dp 3 



;fm(f3m(E-fl)). 



(17) 



(18) 



Starting from the Boltzmann equation we could now derive the equations of motion for the energy density, pressure 
and number density. However since Boltzmann equation is equivalent to the continuity equations ([2| we will only use 
the latter. Therefore by differentiating equation (16 1 with respect to time we get 



C dpp 2 Ef' m fi m (E- M ) = fi n 



C 



fin 

'fin 

fin 
'fir, 

fin 
'fin 

fin 
'fin 



C J™ dp-^(p(E 3 -vE 2 ))f„ 



dpp 

_fin 
fin 



E 2 (E~ f i)df„ 



dp 



dp(E 3 - iiE 2 



3p 2 E-2w 2 )f, 



dp^E(p 2 + to 2 ) - [i{p 2 + to 2 )) f m + 3p m - 2(in r . 



C 



C dpi to (E — n)\ f m + 4p m - 3(in 



(19) 



2/ E-n 



+ 4/9 m - 3(m r , 



where prime now means derivative with respect to the argument of the distribution function and we have also made 
partial integration in the third equality where the assumption /i = m m valid for the dark matter particles was used. 
We can now combine this equation with the continuity equation ^ and solve the equation for fi m . After some algebra 
we get the result 



fin 

fin 



3H( Pm + P rn ) 
Ap m - 3[in m + n 



3ff(l + 

4 - 3 X r, 



) - r g n g /n tr 

m 2 /pr, 



E-jj, 

p 2 



(20) 



i/Pm which is close to zero for relativistic 



Note that in the last equality we have also used a new variable x m = 
matter and one for non-relativistic case. 

The temperature in the early universe is generally very high which means that both Bose-Einstein and Fermi-Dirac 
statistics can be well approximated with the Maxwell-Boltzmann distribution. The distribution of the the dark matter 
particles can be therefore written as i.e. f m = e - ' 5 " 1 ^ - ^. Note that this has also pragmatic reasons that will be 
evident in the numerical calculations. By substituting this distribution in the definition of energy density and after 
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some integral transformations the result can be written in terms modified Bessel functions K n (x): 

= Ce 13 "^ dEE 2 \/ E 2 - m 2 e - f3 ™ E = Ce^m A / dxx 2 ^x 2 - \ e ~^ mx (21) 

Jrn Jl 

2 

= Ce^^fm^K^m^) + 3X 2 (m/3 m )) 

where in the last equality we have used modified Bessel functions of the second kind K 2 (x) . The pressure and number 
density can be written similarly in terms of Bessel functions as 



P m = 



and 



(27r) d p m 



(22) 



(23) 



The equation of state of the dark matter particles is defined as the ratio of pressure to density and by using equations 

K 2 (mf3 m ) 



(21 I and (17) this simplifies to 



L ° m m(3 m Ki{m[3 m ) + 3K 2 (m(3 m ) ' 

We can also write the Xm function in terms of Bessel functions as 

mP m K 2 (mf3 m ) 

m mftinK^mPm) + 3K 2 (m/3 m ) ' 

The final term to be integrated is {{E — fi)/p 2 } which equals now 



E- n 
P 2 



me 



0,1 



(m/3 m K o (m0 m ) + Ki(mf3 m ) 



— m 2 Ki (m0 m ) 
m 



(24) 



(25) 



(26) 



These expressions can now be substituted in the temperature equation ( 20 1 and use some of the identities of Bessel 
functions to get the final result 



(27) 



p m 37J^m/3 m fsri(TO/3 m ) + 4K 2 (mf3 m )j - T a Q< a /Q, m [m[3 m Ki(mf3 m ) + 3K 2 (mf3 m ) 



if3 m (mf3 m K (mf3 m ) + (5 - fif3 m )K 1 (mf3 m )^ + 3(4 - fj,/3 m ) K 2 (m/3 m ) 



which determines the evolution of the temperature of dark matter particles. We have plotted this function in terms of 
Hubble parameter, H, in figure [T] when T m = i.e. when there is no curvaton decay into matter. The figure clearly 
shows that the equations leads to the proper relativistic, T m a a -1 , and non-relativistic, T m oc a~ 2 , limits for the 
temperature evolution. 



We now have a complete set of equations that determine the evolution of the system: namely equations (111, ( 12 1 



(13 1, (24 1, d25k and (27 1 that we can solve for the evolution of the background and the perturbations at first and 



second order. Note that when solving these equations we have used N = ln(a) as the time variable. 



IV. NUMERICAL RESULTS 



We have solved the evolution of the system (11), (12), ( 13 1 , ( 24 1 , ( 25 1 and rt27| with a range of different initial values 



and classified them to physical and non-physical regions i.e. initial values that could have produced the universe we 
observe. We used a similar strategy previously in [43] for a three fluid curvaton decay where the curvaton decays into 
radiation and non-relativistic matter. In order to label different points into physical and non-physical regions we have 
made a range of tests that system needs to pass in order to be classified as physically viable. These are: 
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FIG. 1: Plot of j3 m i '((3 m H) (equation |27l) when there is no curvaton decay into matter as a function of the logarithm of mf3 m - 
When mf3 m C 1 or equivalently T m 3> m $ m /f3 m oc H which corresponds to the usual relativistic limit T m oc a -1 . In the 
non-relativistic case m0 m ^> 1 (or equivalently T m <C m) /3 m //3 m oc 2H which leads to the usual non-relativistic evolution 
T m oc a~ 2 . Equation (27 1 thus clearly leads to physically sound evolution for the temperature. 



1) We check that the curvaton starts to decay before nucleosynthesis. If this condition is not met the universe 
might be radiation dominated at nucleosynthesis but as the curvaton component behaves like a non-relativistic fluid 
it will eventually start to dominate the universe and therefore make the nucleosynthesis non-physical. 

2) The universe should be radiation dominated at nucleosynthesis temperature which has been set to be ~ 10 5 eV 
[5]. In terms of fractional densities this corresponds to [I3J 



n 7 



^ > 10 5 . (28) 

^ cq 



A different choice would naturally lead to different allowed regions but as was noted in [33] the changes are not drastic. 

3) The radiation matter equality should happen at ~ 0.706 eV [42]. This turns out to be quite strict limit in terms 
of the resolution of the simulations we have computed. We have therefore easened this to a range of possible equality 
temperatures: 0.6 eV < T eq < 0.8 eV . This limitation could be circumvented by increasing the number of different 
initial values (and thereby increasing the resolution of the simulation) but this quickly leads to very long simulation 
times. This way we instead get a region which includes the initial values corresponding to the more strict equality 
temperature. 

4) The /nl should be within the limits —9 < Jml < HI given by WMAP 5 year data [3J- The equation for the 
non-gaussianity in the three fluid curvaton decay was derived in |44j and we will only cite the result 

, 25(g 2 -^)-60 gl r 1 +96^-30r 2 

JNL — 7. — 77. i — \ > l^ y j 

6ri(6ri - hqx) 



where 



Scr,m >>a,m ,„_. 

C (2)| A*h 

r 2 = 7^ J 92 = 7Y\ 

V^cr.in/ VscTjin/ 

and where the different numerators are evaluated at the time of decoupling when the un iver se is matter dominated 



and hence £W ~ (j^ . The curvaton initial perturbation values were given in equations (15 1. Note that the system 
is adiabatic if q\ = r\ and qi — r%. The physical meaning of and qi is that they tell how efficiently the initial 
curvature perturbations are converted into the matter and radiation components. 

Before the simulations can be run different initial values have to be set: namely the initial temperatures of radiation 
T x i n and the initial temperature T m - m and particle mass m m for the dark matter, which also define the corresponding 
energy densities, and the initial fractional density of the curvaton. The initial Hubble parameter value can be 
calculated from these variables and the interaction strengths 1^ vary from 10 -30 to 10 5 thereby defining the resolution 
of the simulations. 
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In the usual curvaton hypothesis the curvaton field is usually assumed to be subdominant after inflation and the 
relativistic decay products of inflaton dominate the energy density. The initial temperature of the system is therefore 
closely related to the reheating temperature. The simulations have been run with a different initial temperatures 
corresponding to values T lM = 10 14 , 10 16 , 10 18 and 10 20 eV. 

Since the dark matter particles are produced by the curvaton decay their energy density is assumed to be initially 
small relative to radiation. Note that when starting the simulations equations (20 1 and (27 1 are clearly ill defined 
when Q m — > 0. We have therefore set initially f2 m ^ which can be interpreted that when the simulation starts, 
the decay of the curvaton field has already started. The initial dark matter temperature, T m , has been set to be two 
orders of magnitude smaller in value than the radiation temperature. Other values might lead to numerical problems 
during the evaluations and this selection leads to very small initial dark matter fractional densities, i.e. f2 m .; n <C 1. 

The amount of curvaton in the system is also a free parameter. However since we are following the usual curvaton 
hypothesis initially the curvaton is subdominant. This translates to values CTj i n — 10', i = —4,-3,-2 and — 1 in 
the different simulations. 

Since we have made weak assumptions on the type of the dark matter, the mass m m can have a wide range of 
different values (A review of different dark matter particle models can be found for example in |47) ) . The only 
limitations are that initially dark matter is relativistic i.e. T m 3> m and that it does not interact with radiation. This 
makes it possible to consider a wide range of different dark matter particle species and to study their naturalness in 
the curvaton scenario. 



-15- 




FIG. 2: Contour plot of the system when (a) T lM = 10 18 eV, T mM = 10 16 eV, m m = 10 13 eV and fi CT ,i n = 0.1 and (b) 
Tj,i n = 10 12 eV, T m ,i n = 10 10 eV, m m = 10 9 eV and Q.cr,in = 10 _1 . The plot shows the physicality tests of the system for 
different values of log 10 (T 7 ///o) and \og 10 (T m / Ho). The numbers in different regions indicate the number of passed tests. 
Black points at the edge of white region are the points that passed all the tests, while crosses and circles mark points which 
passed the matter-radiation equality test. Different initial values lead to qualitatively similar results. 



Typical results of the simulations are shown in figures 2(a) and 2(b) The log 10 (r 7 /ifo) and log 10 (r m / H ) plane 
is divided into several separate areas: in the lower left corner the decay widths are too small and the curvaton field 
dominates the evolution of the system even at T eq — 0.7 eV. The triangular area at the lower right corner is caused by 
the non-gaussianity limit. Inside this region r 7 i s much larger than T m and most of the curvaton decays into radiation. 
In terms of non-gaussianity parameter, eq. (29), this means that r± is very small and the generated non-gaussianity 
is very large and therefore labeled as unphysical. The nucleosynthesis limit, eq. (28), is satisfied inside the trapezoid 
shaped region between values log 10 (r 7 /7J ) = — 25 and log 10 (r 7 /i7 ) = —18 in figure 2(a) Points that also satisfy 
the radiation matter equality limit, 0.6 eV < T eq < 0.8 eV , are around the outer edge of this nucleosynthesis region. 
Despite the easened limitations for the equality temperature, this clearly gives only a very narrow region of plausible 
initial values. 

In terms of the naturalness of the curvaton scenario, the decay widths of the curvaton are limited to be very small 
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if the system is compelled to satisfy the previously presented physicality tests. Smaller initial temperatures and dark 
matter particle mass allow larger values but even in this case the decay widths require fine-tuning. 

V. DISCUSSION AND CONCLUSIONS 

We have studied the behavior of a three fluid curvaton decay model when the possibility of more complex thermo- 
dynamics of the cold dark matter decay product have been also taken into account. The assumptions we have made 
are relatively loose: the decay products of the curvaton were assumed to be non-interacting. This is quite a natural, 
if the other fluid is interpreted as conventional radiation and the other as dark matter. The system was simulated for 
a wide range of different initial values including initial temperature, dark matter particle mass and curvaton decay 
widths. This made it possible to study the naturalness of the curvaton model in this model. 

We found that the generated non-gaussianity in the system can be considerable but these large values are reached 
with initial values that do not lead to physically sound evolution for the system. The non-gaussianity is therefore 
expected to be of the same order as in the simple inflaton models [3D] . 

In terms of the dark matter particle mass, m m , the curvaton model leads to physically sound evolution with a range 
of different values but not without fine tuning in the decay widths. Lower initial temperatures and mass m m easen 
this issue but not considerably. We therefore conclude that the current curvaton decay model is not satisfactory from 
the point of view of naturalness. 

According to our analysis, a related, but more viable and physically sound model could be the curvaton decay into 
weakly interacting massive particles (WIMPs) [33]. Interactions between cold dark matter and photons would keep 
them at the same temperature in the early universe until they freeze-out as the temperature drops below their mass. 
The same strategy we used here could be applied also to this scenario to study the naturalness of the curvaton WIMP 
scenario. This is however beyond the scope of this paper and is left as a future work. 
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